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axisymmetric  displacements  are  obtained  for  each  segment  by  solving  a related 
axisymmetric  configuration.  A perturbation  analysis  is  then  performed  to 
match  the  solutions  at  certain  points  between  the  segments,  and  obtain  the 
perturbation  displacements  for  the  total  structure.  The  total  displacement 
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total  displacements  are  known.  The  method  is  applied  to  a number  of  examples 
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I.  INTRODUCTION 


Many  problems  in  structural  analysis  have  an  axisymmetric  configuration. 
However,  there  are  other  problems  which,  while  they  are  almost  axisymmetric, 
have  non-symmetric  features  which  render  purely  axisymmetric  solutions  inap- 
plicable. One  approach  to  solving  these  types  of  problems  would  be  to  use 
a three-dimensional  finite  element  formulation.  However,  if  some  advantage 
can  be  taken  of  the  fact  that  these  geometries  will  yield  solutions  somewhat 
similar  to  the  axisymmetric  cases,  some  efficiency  could  be  gained  over  the 
purely  three-dimensional  approach.  Consequently,  a method  that  makes  use  of 
two-dimensional  axisymmetric  analysis,  but  also  yields  first  order  non-sym- 
metric  effects,  is  an  appealing  alternative. 

This  paper  attempts  to  develop  an  analysis  of  the  stress  and  displace- 
ment field  of  slightly  non-axisymmetric  bodies.  This  is  done  by  obtaining 
the  axisymmetric  displacements  for  various  parts  of  the  body,  and  then 
forcing  the  displacements  to  match  at  certain  points  through  a perturbation 
analysis.  The  axisymmetric  displacements  together  with  the  perturbation 
displacements  combine  to  yield  the  total  displacements  for  each  section  of 
the  body.  Once  the  displacements  are  known,  it  is  a simple  matter  to  cal- 
culate stresses  and  strains  at  any  desired  point  in  the  structure. 
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II.  METHOD  OF  ANALYSIS 


The  purpose  of  the  present  investigation  is  to  develop  an  approximate 
method  of  analysis  of  elastic  configurations  whose  geometry  is  almost 
axisymmetric,  but  has  pronounced  geometrical  changes  in  the  circumferential 
direction.  The  method  of  analysis  which  is  proposed  is  based  on  the  finite 
element  numerical  approach,  and  will  attempt  to  make  some  use  of  the 
geometrical  aspects  of  the  two-dimensional  axisymmetric  analysis . C1) > (2)* 
Although  it  is  not  necessary  for  the  proposed  method,  let  us  for  the  purpose 
of  discussion  assume  that  the  geometry  repeats  at  some  regular  intervals  in 
the  circumferential  direction  and  therefore,  only  one  portion  of  the  total 
structure  need  be  considered.  Of  course,  if  there  is  no  repetitions  then 
the  total  structure  would  be  analyzed.  Consider  now  one  such  portion  as 
illustrated  in  Figure  1.  For  the  purpose  of  analysis  the  configuration 
shown  in  Figure  1 is  divided  into  a set  of  finite  elements  whose  quadrila- 
teral shape  is  defined  in  the  r-z  plane. (2)  At  this  point  direct  use  can  . 

be  made  of  a suitable  mesh  generation  procedure  in  a two-dimensional  plane/ 
On  each  face  of  this  segment  the  finite  element  nodes  are  assumed  to  corres- 
pond to  the  nodes  in  the  adjacent  segments. 

The  first  step  is  to  obtain  the  axisymmetric  solution  for  each  of  these 
segments.  In  addition,  the  stiffness  and  force  matrices  for  the  non-axisym- 
metric  solution  are  calculated  and  stored.  In  general,  the  displacements 
generated  by  the  axisymmetric  solution  will  not  be  compatible  with  the  dis- 
placements generated  for  an  adjacent  segment.  Consequently,  some  changes 
need  to  be  made  in  these  displacements  to  obtain  the  real  situation.  The 
total  displacement  is  now  represented  by  two  components.  If  we  call  the 
displacement  u we  can  write: 

{u}  = {ua}  + {up}  (1) 

where{ua}is  the  displacement  from  the  axisymmetric  solution  and(up)is  the 
additional  perturbation  displacement. 

For  each  side  of  a segment  the  perturbation  displacements  are  expressed 
as  follows: 

ur  = ai  + a2r  + z(aj  + a4r) 

Uz  = a5  + a6r  + z(a7  + a8r)  (2) 

U0  = a9  + a10r  + Z(all  + al2r) 
where  subscripts  r,  z,  0 refer  to  coordinates  in  Figure  1. 
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Superscripts  refer  to  references 


Similar  expressions  can  be  written  for  the  other  side  of  the  segment. 


If  the  total  displacements  are  forced  to  match  for  two  adjacent  segments 
at  four  nodal  points,  a sufficient  number  of  equations  to  be  able  to  solve 
for  the  coefficients  of  the  perturbation  displacement  equations  will  be  ob- 
tained. These  nodes,  where  the  solutions  match  are  called  connecting  nodes. 
Then  for  each  segment: 

{u£  } = [R]  {p}  (3) 

P 

where  (uc  } is  the  perturbation  displacement  at  the  connecting  nodes,  [R]  is 
matrix  of*5 the  connecting  nodal  coordinates,  and  {p } is  the  unknown  perturba- 
tion coefficients.  Solving  for  {p}  yields: 


{p}  = [RJ'1  (uc  } 
P 


Also,  for  each  quadrilateral  element  in  the  segment: 

{up}  = [G]  {p}  (5) 

where  [G ]is  the  matrix  of  nodal  coordinates  for  the  element.  Substituting  equa- 
tion (4)  into  equation  (5)  yields: 

{u  } = [G]  [R]'1  {uc  } (6) 

F P 

Also  for  each  element: 

{fJ  = [Sj]  ({uft}  + {up})  (7) 

Where  [S^Jis  the  non-axisymmetric  stiffness  matrix  and{f}is  the  internal  force 
vector.  Now,  applying  the  principal  of  virtual  work  to  equation  (7)  gives: 

<5UpT  {f } = 6upT  [S1]({ua>  + {up})  (8) 

From  equation  (6),  it  can  be  seen  that: 

T 

T T -1  T 

6u  1 = 6u  [R]  [G]  (9) 

P 

Substituting  this  into  equation  (8)  and  summing  the  contribution  of  all  the 
elements  in  the  segment,  defines  the  following: 


T 

- 1 T 

i rr<i  1 r 


T 
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Wj  = Su  I([R]"X  [G]I[S1]{ua}  + [R]"  [G]'[S  ][R]_1{u  })  (10) 

P P 


Also,  there  is  contributions  from  the  external  forces.  This  can  be  written 


6upT  {fE>  = 6upT  {P1> 


t 


where  {Pj}  is  the  external  force  vector  for  an  individual  element.  The  same 
process  of  substitution  and  summing  over  all  the  elements  defines  the  following 


WE  = 6uc  T l [R]"1  [G]T  {Pl} 
P 


U2) 


Also  there  will  be  a contribution  from  the  concentrated  loads.  This  contribu- 
tion can  be  stated  as: 


V • *U  T {S} 

F P 


(13) 


where  {g}  is  the  concentrated  force  vector  for  the  element.  Then,  for  the 
whole  segment: 


W = <5u  T £ [R]-1  [A]T  {g} 
P 


(14) 


where  [A]  is  a matrix  relating  the  perturbation  displacements  at  the  nodes 
where  the  concentrated  forces  are  applied  to  the  perturbation  coefficient 
matrix.  That  is: 


(up)  = [A]  {p} 

at  nodes  where  concentrated  forces  are  applied. 

Balancing  the  virtual  work  done  by  all  the  forces  gives: 

WE  + WC  = W! 

Here,  the  following  substitutions  can  be  made: 

T 

fe  = [R]"1  (I  [g]t  (Pj)  + [a]t  {g» 


(15) 


(16) 


[SK]  = [R]'1  (Z[G]T[S1][G])[R]_1 


(17) 


{Fj}  = [R]'1  (I[G]T[S1]  {ua}) 

T 

where  the  [R]-l  has  been  brought  outside  the  summation  since  it  is  constant 
for  the  entire  segment,  and  the  summation  is  over  all  the  elements  in  the 
segment.  The  terms  in  the  summation  can  be  easily  calculated  during  the 
axisymmetric  solution  and  later  be  multiplied  by  the  [R ] ~ 1 terms.  Then  {F^}, 

[SK],  and  {Fj}  can  be  stored  for  each  segment.  These  substitutions  give: 


tsK]  {uc  > + (Fj } = {Fp} 

P 


(18) 


6 


Since  it  is  the  total  displacements  that  must  match  at  the  connecting  nodes, 
equation  (18)  can  be  written  as: 


tv  <%>  - [vtuc 1 * fv  ■ 

T a 


(19) 


where  (u 
is  the  c 
except  (u 


} is  the  axisymmetric  displacement  at  the  connecting  nodes  and  {u 
i total  displacement  at  the  connecting  nodes.  Everything  is  known^ 
}.  Combining  terms  in  equation  (19)  yields: 


{ucT}  " {FC} 


(20) 


} 


where  {F^}  is  a combination  of  {F^},  {F^.},  and  [S^]  (ucal.  There  will  be 

one  matrix  equation  for  each  segment.  Now  [S  ] can  be  assembled  into 
a banded  global  stiffness  matrix  for  the  overall^  structure.  Then  appropriate 
boundary  conditions  can  be  applied  to  restrict  the  rigid  body  motion,  preserve 
symmetries,  or  comply  with  external  restraints.  Then  equation  (20)  can  be 
solved  for  the  {u  } vector,  giving  the  total  displacements  at  all  of  the  con- 
necting nodes . CT 


Now  it  is  a simple  matter  to  work  back  from  this  point  to  obtain  the 
stresses  and  strains  for  each  element  in  each  segment.  For  a given  segment 
(u  } and  (u  } are  known  and  therefore,  (u  } can  be  calculated  from: 


{u  } = {u  } - (u  } 
c c„  c 

p T a 

,-l 


(21) 


Then  knowing  {u  } and  [R]  , {p}  can  be  obtained  from  equation  (4) . Then 

for  each  cp  element  the  perturbation  displacements  {u^}  can  be  obtained 

from  equation  (5) 
equation  (1). 


Then  knowing  {u  } and  (u  },  it  is  easy 
a p 


to  get  (u)  from 


At  this  point,  simply  apply  the  base  finite  element  equations: 

(e)  = [B]  {u} 

(22) 

{a}  = [D]  {e} 

where  [B]  and  [D]  were  calculated  and  saved  from  the  axisymmetric  solution, 
to  obtain  the  stresses  and  strains  for  each  element. 


III.  NUMERICAL  RESULTS 


The  first  test  of  this  method  is  that  if  the  body  is  indeed  completely 
axisymmetric,  the  perturbation  analysis  should  give  the  original  axisymmetric 
answers.  This  was  tested  on  two  examples;  a spinning  disk  with  a hole  in  the 
center,  and  a disk  under  internal  pressure.  In  both  cases  the  axisymmetric 
solution  was  duplicated. 

The  next  example  was  a non-axisymmetric  disk  under  internal  pressure. 

The  answers  for  this  example  were  compared  to  answers  obtained  from  a plane 
stress  analysis.  This  example  is  basically  a plane  stress  problem  in  the 


l 


r-0  plane.  The  configuration  that  was  analyzed  is  shown  in  Figure  3.  Due 
to  symmetry  only  one  quarter  of  the  total  body  was  divided  into  six  segments 
for  the  analysis,  numbered  as  shown  in  Figure  3.  The  predominant  stress  for 
this  example  is  cr-g . Results  of  the  perturbation  analysis,  plane-stress 
analysis  and  the  axisymmetric  analysis  for  all  segments  are  shown  in  Table  1. 
It  can  be  seen  that  the  perturbation  analysis  gives  close  agreement  with  the 
plane  stress  solution. 

The  next  example  tested  was  a tube  that  was  axisymmetric  for  half  its 
length  and  had  a cross-section  as  shown  in  Figure  3 for  the  other  half  of 
its  length.  The  only  test  of  accuracy  that  could  be  made  on  this  example 
was  to  see  if  the  answers  tended  to  the  known  limit  solutions  at  each  end. 
That  is  that  near  the  end  of  the  tube  the  answers  for  the  nonaxisymmetric 
part  tended  to  those  presented  in  Table  1,  and  at  the  other  end  the  answers 
tended  toward  the  axisymmetric  solution.  These  results  for  all  segments 
are  presented  in  Table  2.  From  this  we  can  see  that  the  stresses  vary  from 
end  to  end  and  that  they  are  approaching  the  limiting  values  for  the  case 
of  the  tube  having  no  variation  in  the  z direction. 

IV.  CONCLUSION 

From  the  test  examples  used,  it  can  be  concluded  that  this  method  of 
solution  can  handle  non-axisymmetric  geometries.  Also,  the  accuracy  seems 
to  be  sufficient  for  bodies  that  have  geometry  varying  in  the  r-0  plane, 
the  r-z  plane  or  both.  Therefore,  this  method  presents  a simpler  alterna- 
tive to  a purely  three  dimensional  analysis,  and  for  a large  class  of 
problems  this  would  appear  to  have  a great  advantage  in  terms  of  saving  in 
the  computer  execution  time. 
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TABLE  1.  Stress  Ogg  from  various  methods  of  analysis  of  a disk  with 
internal  pressure;  p = 1000  Pa. 


STRESS  or 


Radial  Distance 


6.25 

8.75 

11.25 

13.75 

6.25 

8.75 

11.25 

13.75 

6.25 

8.75 

11.25 

13.75 


Perturbation  Plane  Stress 
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TABLE  2.  Stress  aflfl  from  perturbation  analysis  for  a nonaxisymmetric 
under  internal  pressure  ; p = 1000  Pa. 


STRESS  a00 


Radial 

Distance 

Axisymmetric  end 

Nonaxisymmetric  end 

Disk 

Disk 

6.25 

1 934 

857 

i 

1010 

1164 

8.75 

486 

491 

503 

505 

11.25 

301 

345 

284 

208 

13.75 

205 

271 

108 

46 

6.25 

909 

857 

1043 

1164 

8.75 

484 

491 

505 

11.25 

314 

345 

EH 

13.75 

229 

271 

6.25 

847 

857 

8.75 

480 

491 

511 

505 

11.25 

339 

345 

244 

208 

13.75 

273 

271 

97 

46 

6.25 

799 

857 

1067 

1078 

8.75 

456 

491 

878 

906 

11.25 

328 

345 

13.75 

265 

271 

6.25 

757 

857 

1127 

1078 

8.75 

485 

491 

832 

906 

11.25 

394 

345 

13.75 

354 

271 

6.25 

755 

857 

1155 

1078 

8.75 

479 

491 

811 

906 

11.25 

386 

345 

13.75 

345 

271 

L. 
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APPENDIX  A 


1.  SEGMENT  CONTROL  CARD 

Format  (2110) 

Columns  1-10  NTYPS  (Number  of  different  types  of  segments; 

4 maximum) 

11-20  NTOTS  (Number  of  total  segments;  8 maximum) 

2.  SEGMENT  DATA  CARDS 

One  card  for  each  type  of  segment . 

Format  (F10.5,  110) 

Columns  1-10  THETA  (Angle  subtended  by  segment) 

11-20  NST  (The  number  of  segments  of  each 
type;  5 maximum) 

3.  SEGMENT  NUMBERING  CARDS 

One  card  for  each  type  of  segment . 

Format  (5110) 

Columns  1-10  NUMS  (1)  (Reference  numbers  for  segments  of 

: : each  type  in 

41-50  NUMS  (5)  global  numbering  system;  there  can  be 
one  to  five  numbers  per  card  depending 
on  how  many  segments  there  are  of  each 
type.) 


15 


CONNECTING  NODES  CARDS 

One  card  for  each  segment.  These  cards  must  be  in  order  according  to  the 
global  numbering  system  for  segments. 

Format  (8110) 

Columns  1-10  NPC  (1)  (Nodal  number  for  connecting  nodes 

: : according  to  the  axisymmetric 

71-80  NPC  (8)  grid  for  the  segment) . 

Cards  S-18  must  be  repeated  for  each  different  type  of  segment.  These 
are  the  control  cards  for  the  axisymmetric  solution. 

5.  TITLE  CARD 

Format  (20A4) 

Columns  1-80  TITLE  (Title  for  particular  case) 

6.  CONTROL  CARD 

Format  (615,  F5.0,  515) 

Columns  1-5  NNLA  (Number  of  nonlinear  approximations;  NNLA  = 1 

for  this  version  of  the  program) 

6-10  NUMTC  (Number  of  temperature  cards;  if  -2,  a 
constant  temperature  is  specified) 

11-15  NUMMAT  (Number  of  different  materials;  6 maximum) 

16-20  NUMPC  (Number  of  boundary  pressure  cards;  200  maximum) 

21-25  NUMSC  (Number  of  boundary  shear  cards;  200  maximum) 

26-30  NUMST  (Number  of  boundary  shear  cards  in  tangential 

direction;  200  maximum) 

31-35  TREF  (Reference  temperature) 

36-40  INERT  (This  parameter  decides  if  inertia  loads 

will  be  present,  INERT  = 0 means  zero  values  of  axial 
acceleration,  and  angular  acceleration  and  velocity 
for  each  load  increment) 


1<> 


51-55  INCF  (If  INCF  = 0,  then  surface  loads  for  each  time 

increment  will  be  the  same  as  for  first  increment) 

56-60  IPLOT  (Plot  parameter,  1 if  plot  required) 

7.  MESH  GENERATION  CONTROL  CARD 

Format  (SI5) 

Columns  1-5  MAXI  (Maximum  value  of  I in  mesh;  25  maximum) 

6- 10  MAXJ  (Maximum  value  of  J in  mesh;  100  maximum) 

11-15  NSEG  (Number  of  line  segment  cards) 

16-20  NBC  (Number  of  boundary  condition  cards) 

21-25  NMTL  (Number  of  material  block  cards) 

8.  LINE  SEGMENT  CARDS 

The  order  of  line  segment  cards  is  immaterial  except  when  plots  are 
requested;  in  this  case,  the  line  segment  cards  must  define  the  perimeter  of 
the  solid  continuously.  The  order  of  line  segment  cards  defining  internal 
straight  lines  is  always  irrelevant. 

Format  (3(213,  2F8.3),  IS) 

Columns  1-3  I coordinate  of  1st  point 

4-6  J coordinate  of  1st  point 

7- 14  R coordinate  of  1st  point 

15-22  Z coordinate  of  1st  point 

23-25  I coordinate  of  2nd  point 

26-28  J coordinate  of  2nd  point 

29-36  R coordinate  of  2nd  point 

37-44  Z coordinate  of  2nd  point 

45-47  I coordinate  of  3rd  point 

48-50  J coordinate  of  3rd  point 

51-58  R coordinate  of  3rd  point 

59-66  Z coordinate  of  3rd  point 

67-71  Line  segment  type  parameter 

If  the  number  in  column  71  is 

0 Point  (input  only  1st  point) 

1 straight  line  (input  only  1st  and  2nd  points) 

2 straight  line  as  an  internal  diagonal  (input  only 
1st  and  2nd  points) 

3 circular  arc  specified  by  1st  and  3rd  points  at  the 
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ends  of  the  arc  and  2nd  point  at  the  mid-point 
of  the  arc 

4 circular  arc  specified  by  1st  and  2nd  points  at 
the  ends  of  the  arc  with  the  coordinates  of  the 
center  of  the  arc  given  as  the  3rd  point  (delete 
I and  J for  3rd  point) 

5 straight  line  as  a boundary  diagonal  for  wh  ch  I 
of  1st  point  is  minimum  for  its  row  and/or  I of 
2nd  point  is  minimum  for  its  row  (input  only  1st 
and  2nd  points) 

6 straight  line  as  a boundary  diagonal  for  which  I 
of  1st  point  and/or  2nd  point  is  maximum  for  its 
row  (input  only  1st  and  2nd  points) 

NOTE:  In  specifying  a circular  arc,  the  points  are  ordered  such  that  a counter- 
clockwise direction  about  the  center  is  obtained  upon  moving  along  the 
boundary. 


9.  BOUNDARY  CONDITION  CARDS 

Each  card  assigns  a particular  boundary  condition  to  a block  of  ele- 
ments bounded  by  II,  12,  Jl,  J2.  For  a line  II  = 12  or  J1  = J2.  For  a point 
II  = 12  and  Jl  = J2. 


Format 

(415,  110, 

3F10.0) 

Columns 

1-5 

Mininum  I 

6-10 

Maximum  I 

11-15 

Minimum  J 

16-20 

Maximum  J 

21-30 

Boundary  condition  code 

31-40 

Radial  boundary  condition  code,  XR 

41-50 

Axial  boundary  condition,  XZ 

51-60 

Tangential  boundary  condition  XT 

I f the 

number  in 

Columns  21-30  is 

XR  is 

the  specified  R-load  and 

0 

XZ  is 

the  specified  Z-load  and 

XT  is  the  specified  T-load 

XR  is  the  specified  R-displacement  and 
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i 

XZ 

is 

the 

specified 

Z-load  and 

XT 

is 

the 

specified 

T- load 

XR 

is 

the 

specified 

R-load  and 

2 

XZ 

is 

the 

specified 

Z-displacement 

and 

XT 

is 

the 

specified 

T- load 

1 

XR 

is 

the 

specified 

R-displacement 

and 

3 

XZ 

is 

the 

specified 

Z-di splacement 

and 

XT 

is 

the 

specified 

T-load 

4 

XZ 

is 

the 

specified 

Z-load  and 

XT 

is 

the 

specified 

T-displacement 

XR 

is 

the 

specified 

R-displacement 

and 

5 

XZ 

is 

the 

specified 

Z-load  and 

XT 

is 

the 

specified 

T-displacement 

XR 

is 

the 

specified 

R-load  and 

1 

6 

XZ 

is 

the 

specified 

Z-displacement 

and 

XT 

is 

the 

specified 

T-displacement 

s 

XR 

is 

the 

specified 

R-displacement 

and 

! 

7 

XZ 

is 

the 

specified 

Z-displacement 

and 

K 

XT 

is 

the 

specified 

T-displacement 

NOTE:  All  loads  are  considered  to  be  total  forces  acting  on  one  radian  segment. 


10.  MATERIAL  BLOCK  ASSIGNMENT  CARD 

Each  card  assigns  a material  definition  number  to  a block  of  elements 
defined  by  the  I,  J coordinates. 

Format  (515,  2F10.0,  215) 

Columns 


1- 

-5 

Material 

definition 

number 

(1  through  6) 

6- 

-10 

Minimum 

1 

11- 

-15 

Maximum 

I 

16- 

-20 

Minimum 

J 

21- 

-25 

Maximum 

J 

26- 

-35 

Material 

principal 

property 

inclination 

angle 

BETA 

in  R-z  plane 

36- 

-45 

Materia  1 

pr i nc i pa  1 

p ropertv 

inclination 

angle 

ALPHA 

in  N-T  p 

lane 

10 
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46-50  IANG  (If  IANG  = 0,  then  ALPHA  is  same  for  total 

material  block.  If  IANG  = 1,  the  ALPHA  varies 
in  sign  in  the  I direction  from  element  to 
element  every  NANG  elements.  This  will  allow 
for  equal  but  opposite  helical  angles.) 

51-55  NANG  (Number  of  elements  in  the  I direction  with 
the  same  ALPHA) . 


11.  PLOT  TITLE  CARD* 


Format  (20A4) 
Columns  1- 


Title  (Title  printed  under  each  plot) 


12.  PLOT  GENERATION  INFORMATION  CARD* 

Format  (2F10.0) 

Columns  1-10  RMAX  (Maximum  r coordinate  of  mesh) 
11-20  ZMAX  (Maximum  z coordinate  of  mesh) 

*N0TE : Use  only  if  IPLOT  = 1 (plot  required) 


13.  TEMPERATURE  FIELD  INFORMATION  CARDS 

If  NUMTC  in  columns  6-10  of  the  CONTROL  CARD  is  greater  than  1, 
the  temperature  field  is  given  on  cards.  One  card  must  be  supplied  for  each 
point  for  which  a temperature  is  specified. 

Format  (3F10.0) 


Columns 


11-20 


R coordinate 
Z coordinate 


21-30  Temperature 

If  NUMTC  in  columns  6-10  of  the  CONTROL  CARD  is  -2,  a constant  temperature  field 
is  specified;  the  value  is  given  on  a single  card. 

Format  (F10.0) 

Columns  1-10  Temperature 
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14.  MATERIAL  PROPERTY  INFORMATION  CARDS 

The  following  group  of  cards  must  be  specified  for  each  material 
(maximum  of  6) . 

a.  MATERIAL  IDENTIFICATION  CARD 

Format  (215,  2F10.0) 

Columns  1-5  Material  identification  number 

6-10  Number  of  temperatures  for  which  properties 
are  given  (12  maximum) 

11-20  Mass  density  of  material  (if  required) 

21-30  Thermal  expansion  parameter  (If  1,  free 

thermal  expansions  on  the  material  property 
cards;  otherwise,  coefficients  of  thermal 
expansion  are  on  the  material  property  cards.) 

b.  MATERIAL  PROPERTY  CARDS 

Two  cards  are  required  for  each  temperature. 

First  Card 
Format  (7F10.0) 

Columns  1-10  Temperature 

11-20  Modulus  of  elasticity,  E., 

21-30  Modulus  of  elasticity, 

31-40  Modulus  of  elasticity,  Eq 

41-50  Poisson's  ratio,  v.,„ 

No 

51-60  Poisson's  ratio,  v ^ 

61-70  Poisson's  ratio,  v„Q 

bu 

Second  Card 
Format  (6F10.0) 

Columns  1-10  Shear  Modulus, 

11-20  Shear  Modulus, 

21-30  Shear  Modulus,  G~„, 

oN 

31-40  a T or  a 


n n 


15.  INERTIA  LOAD  CARD 


Format  (3F10.0) 

Starting  with  this  input  card  and  including  the  boundary  force  cards, 
this  data  is  to  be  inputted  as  a block  for  each  load  step,  that  is 
NLINC  times.  There  are  the  following  exceptions  to  this: 


a)  If  INERT  = 0,  then  this  card  is  to  be  omitted  completely 
(no  inertia  load) . 

b)  If  INCI  = 0,  then  this  card  is  not  repeated,  but  appears 
in  first  block  only  (the  inertia  loads  are  constant  for 
each  load  step) . 

c)  If  INCF  = 0,  then  the  following  boundary  pressure  and  shear 
cards  are  to  be  given  only  for  the  first  block  and  not  re- 
peated again  (the  pressure  and  shear  loads  are  constant  for 
each  load  increment) . 

Columns  1-10  ACELZ  (axial  acceleration 

11-20  ANGVEL  (angular  velocity) 

21-30  ANGACC  (angular  acceleration) 


16.  BOUNDARY  PRESSURE  CARDS 


One  card  is  required  for  each  boundary  element  which  is  subjected  to 
a normal  pressure,  that  is  the  number  of  these  cards  is  NUMPC  for 
each  load  increment. 

Format  (215,  F10.0) 

Columns  1-5  Nodal  point  M 
6-10  Nodal  point  N 

11-20  Normal  pressure 

As  shown  in  the  figure  below,  the  boundary  element  must  be  on  the  left 
when  progressing  from  M to  N.  Surface  normal  tension  is  input  as  a 
negative  pressure. 


17.  BOUNDARY  SHEAR  CARDS 


One  card  is  required  for  each  boundary  element  which  is  subjected 
to  surface  shear,  that  is,  the  number  of  these  cards  is  NUMSC  for 
each  load  increment. 

Format  (215,  F10.0J 
Columns  1-5  Nodal  point  M 
6-10  Nodal  point  N 

11-20  Surface  shear 

As  shown  in  the  figure  below,  the  boundary  element  must  be  on  the 
left  when  progressing  from  M to  N.  The  positive  sense  of  the 
shear  is  from  M to  N. 


18.  BOUNDARY  TRANSVERSE  SHEAR  CARDS 


One  card  is  required  for  each  boundary  element  which  is  subject 
to  transverse  shear,  that  is,  the  number  of  these  cards  is 
NUMSC  for  each  load  increment. 

Format  (215,  F10.0) 

Columns  1-5  Nodal  point  M 
6-10  Nodal  point  N 

11-20  Surface  transverse  shear 


M 


19.  BOUNDARY  CONDITION  CONTROL  CARD 


Format  (15) 

Columns  1-5  NRDF  (Number  of  boundary  conditions) 


20.  BOUNDARY  CONDITION  CARDS 

There  are  NRDF  of  these  cards . 

Format  (110,  F10.0) 

Columns  1-10  NREQ  (The  number  of  the  equation  to  be  modified 

in  the  assembled  matrix) . 

11-20  U (The  actual  boundary  condition  valve). 
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APPENDIX  B 


File  Number 
1-2 
3 


7-14 


17 


21-24 
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Useage 

Used  in  axisymmetric  solution 

T 

Saves  product  of  [G]  [S^]  for  each  element 

from  one  part  of  axisymmetric  solution  for  later 
use. 

These  files  are  used  sequentially  for  each 
different  segment  (i.e.  data  for  segment  1 is 
stored  in  file  7;  for  segment  2 in  file  8,  etc.). 
Data  saved  are  { FE } , {Fl},  [SK]  and  {u^ } , for 
use  in  the  perturbation  analysis 

Saves  the  axisymmetric  displacement  vector  and 
nodal  point  connections  for  each  type  of  segment 
(i.e.  data  for  segments  type  1 is  stored  in  file 
17  first.  Then  data  for  segment  type  2 is  stored 
in  file  17,  etc.). 

Saves  [G],  [CRZ] , and  [BS^]  for  each  element  of 
the  segment.  Data  for  each  different  type  of 
segment  goes  in  a different  file.  This  data 
is  used  to  calculate  stresses  and  strains  when 
the  total  displacements  are  known. 

Saves  [R^ ] 1 for  each  different  segment  in  order. 
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